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Abstract We provide strong Lp-rates of approximation of nonsmooth integral-type function¬ 
als of Markov processes by integral sums. Our approach is, in a sense, process insensitive and 
is based on a modification of some well-developed estimates from the theory of continuous 
additive functionals of Markov processes. 
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1 Introduction 

Let Xt, f > 0, be an R'^-valued Markov process. We study an integral functional 

^T{h) = [ h{Xt)dt 

JQ 

of this process. The most natural numerical scheme to approximate such a functional 
is the sequence of integral sums 

rj^ n 1 

dT,n{h) = — ^ h{X(kT)/n), n>l, 

^ k^O 
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and the main objective of this paper is to study approximation rates within this scheme. 
The function h, in general, is not assumed to be smooth, and therefore the mapping 



may fail to be Lipschitz continuous (and even simply continuous) on a natural func¬ 
tional space of the trajectories of X (e.g., C{0,T) or D{0,T)). This makes it im¬ 
possible to carry out the error analysis with a classical technique (see, e.g., [5]). The 
typical case of interest here is h = 1 a, with /t(^) being respectively the occupation 
time of X at the set A up to the time moment T. 

In the paper, we establish strong Lp-approximation rates, that is, the bounds for 


E\lT{h) - lTAh)\'’ 


Our research is strongly motivated by the recent paper [7], where such a problem was 
studied in a particularly important case where X is a one-dimensional diffusion, and 
we refer the reader to [7] for more motivation and background on the subject. The 
technique developed in [7], involving both the Malliavin calculus tools and the Gaus¬ 
sian bounds for the transition probability density, relies substantially on the structure 
of the process, and hence it seems not easy to extend this approach to other classes 
of processes, for example, multidimensional diffusions or solutions to Levy driven 


SDEs. 


We would like to explain in this note that, in order to get the required approxima¬ 
tion rates, one can modify some well-developed estimates from the theory of continu¬ 
ous additive functionals of Markov processes. An advantage of such approach is that 
the assumptions on the process are formulated only in terms of its transition proba¬ 
bility density and therefore are quite flexible. The basis for the approach is given by 
the fact that the weak approximation rates for 


Elrih) - ElT,n{h) 


are available as a consequence of a bound for the derivative w.r.t. t of the transition 
probability density; see [4], Theorem 2.5, and Proposition 2.1 below. To explain the 
principal idea of the approach, let us assume for a while that h is nonnegative and 
bounded. Then the integral functional Irih) is a W-functional of the process X\ 
see [1], Chapter 6. It is well known that the properties of a PL-functional are mainly 
controlled by its characteristic, that is, the expectation 


Exlrih). 


In particular, the convergence of characteristics implies the L 2 -convergence of the 
respective functionals. The core of our approach is that we extend the Dynkin’s tech¬ 
nique for a study of convergence of PL-functionals and give approximation rates 
for integral functionals Irih) by difference functionals lT,n{h), based on the weak 
approximation rates for their expectations. We remark that now we are beyond the 
scopes of the original Dynkin’s theory because Irih) may fail to be a PL-functional 
(we do not assume h to be nonnegative), and lT,n{h) definitely fails to be a PL- 
functional. In addition, Dynkin’s theory addresses L 2 -bounds, whereas, in general. 
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we are interested in Lp-bounds. This brings some extra difficulties, which however 
are not really substantial, and we resolve them in a way similar to the one used in the 
classical Khas’minskii lemma; see, for example. Lemma 2.1 in [8]. 

2 Main results 

2.7 Notation, assumptions, and auxiliaries 

In what follows, Pj; denotes the law of the Markov process with Xq = x, and de¬ 
notes the expectation w.r.t. this law. The natural filtration of the process X is denoted 
by {Xt,t > 0}. The process X is assumed to possess a transition probability density, 
denoted below by pt{x, y). By C we denote a generic constant; the value of C may 
vary from place to place. Both the absolute value of a real number and the Euclidean 
norm in are denoted by | • |. 

Our standing assumption on the process X under investigation is the following. 
X. The transition probability density pt{x, y) is differentiable w.r.t. t and satishes 


Mx,y) < - y)), t<T, 

dtpt{x,y)\ < - y)), t < T, 


( 1 ) 

( 2 ) 


with some fixed a G (0, 2] and distribution density Q. 

The assumption X is motivated by the following class of processes of particular 
interest. 

Example 2.1. Let XT be a symmetric a-stable process with a G (0, 2]; in the case 
a = 2 this is just a Brownian motion. Then 




with being the distribution density of Xi. Respectively, (2) holds with Q = Qa, 



where C 2 < {2EXi)~^ and ci, c should be chosen such that Q{x) dx = 1. 

Observe that, in a sense, this bound is “stable under perturbations of the process 
X.” Namely, if X is a uniformly elliptic diffusion with Holder continuous coeffi¬ 
cients, then (2) with Q = Q 2 and properly chosen C 2 is provided by the classical 
parametrix method', see [3]. An analogue of the parametrix method for a-stable gen¬ 
erators with state-dependent coefficients yields the bound (2) with Q = Qa', a' < a, 
for a-stable driven processes X; see [2] and [6]. 

Our principal assumption on the function h is the following. 

HI. The function h satisfies 


where V : R'*' —>• [1, -foo) is a hxed function such that 
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• fga V(T\x\)Q{x) dx < oo, T > 0; 

• V{ri) < V{r2), ri < r2\ 

• y is submultiplicative, that is, 

V{ri + r2) <V{ri)V{r2), ri,r2eK+. 


Observe that for a bounded h, condition HI holds trivially with V = 1. On the other 
hand, in particular cases, one can weaken the assumptions on h by using nontrivial 
“weight functions” V. For instance, if Q = Qoc from the above example, then one 
can take 

V{r) = I^ f ^ R+, 

\(l + r)^ as (0,2), 

with arbitrary C and /3 £ (0, a). We denote 

ii.ii \hix)\ 

The following auxiliary statement is crucial for the whole approach. Its proof is com¬ 
pletely analogous to the proof of (a part of) Theorem 2.5 [4], but in order to make the 
exposition self-sufficient, we give it here. 

Proposition 2.1. Let X and HI hold. Then 

\ExlT{h) — Ea:lT.n{.h)\ < ^ ^ ^^T,Q,v\\h\\vV {\x\) 

Ct,q,v = Tmax|cT(^^_^ Q{y)V{T^/^\y\) dy 


with 


1 . 


Proof. Write 


rTjn f. rp 

ExlT{h) - ExlT,n(h) = / / pt{x,y)h{y)dydt - h{x) 

Jo ^ 

{pt{x, y) - P{k-i)T/n{x, y))h{y) dydt. 


fkTIr. 


^_2 J — jn ^ 

We have, by the bound forp((a;, y) in (1) and properties of V, 


fT/r. 


/R'* 


Pt{x,y)h{y) dydt 


[Tin 

< Crllhlly / t 
Jo 

T j 

= CtWHv [ 


—djc 




Q{t ^/°‘{x-y))V(\y\)dydt 


/R'i 


Q{z)V(\x + dzdt 


<n-^TCT\\h\\vV{\x\) / Qiz)V{T^^‘^\z\)dz. 


( 3 ) 
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Next, 


fkTjn 


k=2 '^{k-i)T/n . 


{pt{x,y) - P(k-i)T/n{x,y))h{y)dydt = / Kn,T{x,y)h{y) dy 


with 


't- pkTjn p 

Kn,T{x,y) ='f 2 / 

J (k—DT/n J 0 


Then 


k=2 ik—i)T/n J (k—l)T/n 
n 

k=2 ' 

T I 

\Kn,T{x,y)\ < — / \dsPsix,y)\ds, 

n Jt/u 


dsPs{x,y) dsdt 


rkT/n 


J{k-l)T/n 

V n J 


Tjn 

and therefore, using the bound for dtPt{x, y) in (2), we obtain, similarly to (3), 


Kn,T{x,y)h{y) dy 


< -CtWHyT s-^[ s-‘^/^Q{s-^^^{x-y))V{\y\)dyds 

^ JT/n Jr<‘ 

< ^CT\\h\\vV{\x\) (^I^Q{y)V{T^/^\y\) dy'^ ds 


= -{logn)CT\\h\\vV{\x\) ( Q{y)V{T^/^\y\) dy ), 


which completes the proof. 


□ 


2.2 Approximation rate in terms o/1|/i||y 

Our main estimate, in a shortest and most transparent form, is presented in the fol¬ 
lowing theorem, which concerns the case where the only assumption on h is that the 
weighted sup-norm \\h\\v is finite. 

Theorem 2.1. Let X and HI hold. Then for every p > 2 such that 




VP{T^^°‘\x\)Q{x) dx < CX), 


we have 


EMT{h)-lTAh)f <C 


logn\ 


p/2 


rvvp{\x\) 


with constant C depending on T, Q, V,p only. 

Proof. Denote, for t € [kTjn, {k + l)T/n), 

/■ u\ + 1)^ 

Pn(t) = -, Qn[t) = -, 

n n 

and write the difference It{h) — It,n{h) in the integral form: 


■= h{h) - It,n{h) = [ Anis)ds, (s) :=/i(Xs) -(s)). 

^0 
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Hence, this difference is an absolutely continuous function of t, and using the Newton- 
Leibnitz formula twice, we get 

\jT,n{h)\^ =p(j)-l) J \js,n{h)\^ An{t)dt^ds. 


We then write 

IP 


I 1^ = Pb- 1) ,n,p(^)) < Pb” 1) (^T.n,p(^) +-ffT.n.p W) , 

where 

HT^n,pih)= \js,n{h)f~^A^is)(^J^’^^\nit)dtys, 

fT rCn{s) 

HT,n,p{h)= \Js,n{h)f \A^{s)\ A^{t)dtds, 

Jo Js 

H^npih)= [ \Js,u{h)r^^n{s)( [ A^{t)dt]ds. 

Jo \J(n.{s) J 

Let us estimate separately the expectations of H^^p{h) and By the 

Holder inequality, 

E^Hy^Jh) < \J,,r^{h)fds 

Again by the Holder inequality, 

I A /„MP/2 


1-2/p 


rCn(s) 


A„(t) dt 


p/2 ^ 2/p 

ds 


E^r. 


KAs) 


A„(f) dt 


p/2 


< 


p/2-1 „T Kr.{s) 


ds 

I p/2 




0 Js 


Because t G [s, Cn(-5)], we have rin{t) = and, consequently. 


£;,|A„b)r^"|A„(f)| 

= E,\h{X,) - 

< {E,\h{X,) - h{X^^(^s))\Y\E,\h{X0 - 

By the properties of V and the bound (1) we have 

E,VP{\Xr\) = [ Pr{x,y)VP{\y\)dy 


I p/2 


( 4 ) 








Rates of approximation of non-smooth integral type functionals of Markov processes 


123 


<CTr-^/^ [ Qir-^^^{x-y))VP{\y\)dy 

< CtVP{\x\) Qiy)VP{T^/‘^\y\) dy^ 

for any r S (0, T], Using this bound with r = t,s, r]n{s) and recalling that |Cn(s) — 
s\ < 1/n, we get 

(^^y\h\\lvmx\) 

with constant C depending on T, Q, y,p only. 

Next, observe that, for every s, the variables 

A„(S), \Js,nih)r^^n{s) 


are -measurable. Hence, 

ExH^^n.pW = Ex(^J \Js,n{h)\^ ‘^An{s)E^^J A„(f) df | (s)^ ds 

<E^( [ |Js,„(d)|^ ^|A„(s)| E^f [ An{t) dt\E^^(^s) 

\do Vdc„(s) 


By Proposition 2.1 and the Markov property of X we have 


ds I. 




< C 


/ logi 

\ n 


k^(l^Cn(s)l)- 


Hence, again, using the Holder inequality we get 

,/logn'' ^ 


Similarly to (4), we have 

£;,|a„(s)|^/Vp/2(|x^„(,)|) < cvp{\x\)\\hr^\ 


1-2/p 

2/p 


Hence, the above bounds for E^H^ „ p{h) and E^H^ ^ ^{h) finally yield 

Ex\JT,n{h)\'’ <c[e, |j,.„(d)|'’ds^ (^^y\h\\lv'^{\x\) (5) 

with a constant C depending on T,Q,V,p only. It can be seen easily that in this 
inequality one can write arbitrary t < T instead of T, with the same constant C. 
Taking the integral over t, we get 


E^ 


\Jt,u{h)\'’ dt < Ct(E; 


\JsAh)fds 



0 
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Because \\h\\v < oo and satisfies the integrability condition from the condition 
of the theorem, the left-hand side expression in the last inequality is hnite. Hence, 
resolving this inequality, we get 

\Js,n{h)fds<{CT)P/^(^-^Y \\hfyVP{\x\), 

which, together with (5), gives the required statement. □ 


2.3 An improved approximation rate for a Holder continuous h 

In this section, we consider the case where h has the following additional regularity 

property. 


H2. The function is Holder continuous with index 7 G (0,1], that is. 


II/1II7 := sup 

x^V 


\h[x) - h{y)\ 

\x-yV 


< 00 . 


An additional regularity of h allows one to improve the accuracy of the previous 
estimates. Namely, the following statement holds. 

Theorem 2.2. Assume thatX, HI, andHl hold. Then, for every p > 2 such that 

f \xpPQ{x) dx < 00 

and 

f vP/'^{TE‘^\x\)Q{x)dx < 00, 

Jr<‘ 

we have 

E,\lT{h)-lT,n{h)\^ < ri-('^P)/(2“)||h||P/2(||h||P/2+||/,||P/VP/2(|x|)) 

with constant C depending on T, Q, V,p, 7 only. 

Proof. The method of the proof remains the same as that of Theorem 2.1 ; hence, we 
use the same notation. The only new point is that, instead of the bound 

now a more precise inequality is available, based on the Holder continuity of h. 
Namely, we have 

= E^\h{x,) - 

< ||h||PT;,|X,-X^„(«)|^'’<C||h||P|s-77„(s)|'^^/“, (6) 

where C depends on T, Q, p, 7 only. 

The last inequality holds due to the following representation. By the Markov prop¬ 
erty of X, for r < s, we have 


E,\X,-Xr\^P = EJ{Xr) 
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where 

f{z)=[ ps-r{z,y)\y-zl'^Pdy 

<Ct [ {s-r)-‘‘/^Q{{s-r)-^/^{z-y))\z-yrPdy 

= CT{s-ryP/‘^ [ Q{y)\yrPdy. 

Jr<‘ 

Thus, for f G [s, CnCs)]^ we have 
= E,\h{Xs) - 

< {E,\h{Xy - h{X^^^s))\^ffE,\h{Xt) - 

and 

/ fT \ 1-2/p / , X l+27/a 

E,Hy^^yh)<ci^E,J^ \JsAh)\^dsj (^-j ||/i||2 

with constant C depending on T, (5,p, 7 only. 

Next, using ( 6 ) and (4), we have 

S,|A„(s)|p/Vp/"(|A^„(,)|) < C\\h\\P/^n-^^P'>/^^^WP/^\x\) 

and 

EooHl„Jh)<c(^E, \Js,n{h)fds^ (^-2^'^n-'^/^\\h\\v\\h\\^V{\x\). 

Hence, the previous bounds for E^H^ „ p{h) and E^H^ npW finally yield 
Ex\jT,n{hf 

<c(^E, £\Js,n{h)fds'^ 

(7) 

with constant C depending on T, Q, V,p, 7 only. Using the same procedure as that 
provided at the end of the proof of Theorem 2.1, we simply have 



which, together with (7), completes the proof. □ 
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2.4 Discussion 

The results of Theorem 2.1 and Theorem 2.2 should be compared with Theorem 2.3 
in [7], where, in our notation, the following bounds were obtained; 


E,\lT{h)-lTAh)\'’< 


C'„-p(1+7)/2^ 

C'(logn)Pn-P(i+T')/2, 


7e (0,1), 

7 = 1. 


In Theorem 2.3 of [7], h satisfies our conditions HI and H2 with 17(|a;|) = 
and X is a one-dimensional diffusion with coefficients that satisfy some smoothness 
condition (assumption (H)); in particular, X holds with a = 2. In this case, our bound 

C'(logn)Pn-P/2-(P7)/4^ 


given in Theorem 2.2, is somewhat worse. On the other hand, this bound is of an 
independent interest because of a wider class of processes X it applies to. 
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